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ABSTRACT 

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength 
of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program 
and its role in the complement of Herschel surveys, and the field selection which includes popular multiwavelength fields such as 
GOODS, COSMOS, Lockman Hole, ECDFS, EGS. We provide an account of the observing strategies and data reduction methods 
used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high 
redshift galaxy populations, thus testing and superseeding previous extrapolations from other wavelengths, and enabling a wide range 
of galaxy evolution studies. 

Key words. Surveys - Galaxies: evolution - Galaxies: active - Infrared: galaxies 
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average can be constrained by comparing the observed ultravi- 
olet emission to the sum of ultraviolet and infrared emission 
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Table 1. PEP fields 
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Notes. List of the six blank deep fields, ten lensing clusters and two z~l clusters to be observed by PEP. As of early June 2011, surveys for all 
fields except Abell 370 and Abell 1689 are completed (see column 'Observations'). Field sizes are quoted in terms of the nominal scan length. 
Given the 3.5'xl.75' size of the PACS arrays, coverage is unequal at the field edges starting to decrease already inside the quoted fields, but also 
extending beyond them. Observing times are as computed in HSPOT 5.2, true execution times can differ due to overhead changes during the 
Herschel mission. 

{a) Approximate. See Sect.|3]for detailed implementation. 

<b) 250, 350, 500 /urn obtained in coordinated observations by the HerMES key program (Oliver et al. 1201 lb . 
(c) In array coordinates, true position angle on sky will depend on execution date of the observation. 



instruments, this has changed dramatically. With its 3.5 m pas- 
sively cooled mirror it provides the much improved spatial reso- 
lution (thus reduced source confusion) and the sensitivity needed 
for the next significant step in far-infrared studies of galaxy 
evolution. Members of the PACS instrument consortium, the 
Herschel Science Centre and mission scientist M. Harwit have 
joined forces in the PACS Evolutionary Probe (PEP) deep ex- 
tragalactic survey, to make use of this opportunity. PEP aims to 
resolve the cosmic infrared background and determine the na- 
ture of its constituents, determine the cosmic evolution of dusty 
star formation and of the infrared luminosity function, elucidate 
the relation of far-infrared emission and environment and deter- 
mine clustering properties. Other main goals include study of 
AGN/host coevolution, and determination of the infrared emis- 
sion and energetics of known high redshift galaxy populations. 

PEP encompasses deep observations of blank fields and lens- 
ing clusters, close to the Herschel confusion limit, in order to 
probe down to representative high redshift galaxies, rather than 
being restricted to individually interesting extremely luminous 
cases. PEP is focussed on PACS 70, 100, and 160 fim observa- 
tions. SPIRE observations of the PEP fields are obtained in co- 
ordination with PEP by the HerMES survey (Oliver et al. 1201 11 1. 
Larger and shallower fields are observed by HerMES (70 deg 2 ) 
as well as by the H- ATLAS survey (570 deg 2 , Eales et al.l20TOt . 
while the GOODS-Herschel program (Elbaz et al. 1201 It pro- 
vides deeper observation in (part of) the GOODS fields that 
are also covered by PEP. Finally, the Herschel lensing survey 
(Egami et al. 2010) substantially increases the number of lens- 
ing clusters observed with Herschel, adding about 40 clusters to 
the 10 objects covered by PEP. Fig.[T]compares for 160/im wave- 
length the area and exposure of the PEP surveys (Table [TJ with 
that of these other major Herschel extragalactic surveys. 



In this paper, we describe the field selection, observing strat- 
egy and data analysis methods of PEP. We give a complete 
overview of the planned PEP observations and their execution 
status as of June 2011 (Table[T]l. We provide a detailed account 
of the science demonstration phase (SDP) data sets for GOODS- 
N and Abell 2218, and give an overview of first science results. 

2. Field selection 

A key element in selecting a field for a deep Herschel extra- 
galactic survey is the availability of a strong multi-wavelength 
database from X-rays to radio wavelengths, which is fundamen- 
tal to many of the science results discussed in Sect. [6] In particu- 
lar, deep optical, near-IR, and Spitzer imaging, as well as a com- 
prehensive set of photometric and spectroscopic redshifts are an 
essential asset for most of the far-infrared studies of galaxy evo- 
lution that we envisage. Deep X-ray data are invaluable for using 
the potential of Herschel for studying the AGN - host galaxy co- 
evolution. 

Another requirement is a low galactic far-infrared back- 
ground in order to minimize contamination by galactic 'cirrus' 
structure and by individual galactic foreground objects. This 
naturally coincides with the selection criteria of extragalactic 
surveys at other wavelengths. In the X-ray regime, for exam- 
ple, these are pushing for a low galactic foreground obscura- 
tion. Given that the power spectra of far-infrared emission of 
galactic cirrus are steeply decreasing towards higher spatial fre- 
quencies (e.g., Kiss et al. 120011 1. such cirrus contamination is 
much less of a practical worry for point source detection in 
deep Herschel fields compared to the previous smaller cryo- 
genic space telescopes. In any case, our blank field selection in- 
cludes some of the lowest galactic foreground fields (Lockman 
Hole, CDFS, 100 /mi sky brightness ~0.4 MJy/sr), reaching up 
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Fig. 1. Area and 160/vm exposure of PEP surveys (red) com- 
pared to other major Herschel extragalactic surveys. Exposure 
is defined as survey time (including overheads) multiplied with 
the ratio of PACS array size and survey area. 1000 sec exposure 
corresponds to an approximate 3cr depth of 8 mJy. Scaling this 
as exposure -0 5 to smaller/deeper fields will be too optimistic 
because of increasing overhead fractions and beginning source 
confusion. 



to -0.9 MJy/sr (COSMOS). Our individual cluster fields typi- 
cally show a 100 fim galactic foreground of 1-2 MJy/sr, with a 
maximum of 4.5 MJy/sr (Abell 2390). 

Our largest field is the contiguous 2 square degree COSMOS 
(Scoville et al. 2007 ) observed for about 200 hours to a 3<x depth 
at 160 fim of 10.2 mJy. At this level, integral number counts 
reach one source per 24 beams (Berta et al. 20101 [20TTb . similar 
to the 5<x 40 beams/source definition of the 'confusion limit' 
used by, e.g., Rowan-Robinson et al. (2001). Slightly deeper 
observations have been obtained for the Lockman Hole (e.g. 
Hasinger et al. 2001), the Extended Groth Strip EGS covered by 
the Aegis survey (Davis et al. 2007), and the extended Chandra 
Deep Field South (ECDFS) for which the name has been coined 
in the X-rays (Lehmer et al. 2005) but a multitude of data ex- 
ist at other wavelengths. Refined analysis of source confusion 
(Dole et al. I20041 > suggests that with the 100 fim and 160 fim 
number counts turning over at a depth that is reached with PEP 
(Berta et al. 120101 1201 11 1, deeper observations can still be ex- 
tracted reliably, in particular when using position priors from 
very deep 24 fim data (e.g. Magnelli et al. 2009). We make use 
of this in our deepest blank field observations which are cen- 
tered on the GOODS fields (Dickinson et al. in prep.), with 
strong emphasis on the GOODS-S. GOODS-S also is the only 
field which we additionally observe at 70 fim, to a depth where 
Spitzer/MlPS would be confused. For redshifts around 1, our 
blank field observations sample a range of environments from 
the field to moderately massive clusters, which are known in par- 
ticular in COSMOS due to its large angular size but still excel- 
lent multiwavelength characterisation. In order to extend this to 
the full range of environments at this redshift, we add dedicated 



observations of two of the best studied massive z~l clusters, 
RXJ0152.7-13 57 (e.g . Ebeling et al. [2000T i and MS1054.4-0321 
(e.g.Tranetal. [T999l . 

We use the amplification provided by massive galaxy clus- 
ters to study sources that would otherwise be too faint for di- 
rect observations even in our deepest blank fields, and to provide 
highest quality SEDs on more luminous objects that are also de- 
tectable in the blank fields. For that purpose, we have selected 
ten of the best studied lensing clusters at redshifts z~0.2-0.5. 
With 4' size our observations for these fields are optimized for 
studying lensed background objects. They will additionally de- 
tect part of the cluster members from the central region, but not 
the cluster outskirts and the adjacent field population. 

Main parameters of our survey fields are summarized in 
TableQ] 

3. Observing strategy 

The scan map is the PACS photometer observing mode that is 
best suited to get deep maps close to the confusion limit, for 
regions that are dominated by (almost) point sources. During a 
PACS prime mode scan map observation, the telescope moves 
back and forth in a pattern of parallel scan lines that are con- 
nected by short turnaround loops. During such a scan, the PACS 
photometer arrays take data samples at 10 Hz frequency. Given 
the presence of 1/f noise in the PACS bolometers, we consis- 
tently follow the recommendation to adopt a medium scan speed 
of 20 arcsec/sec that reduces the effects of 1/f noise on point 
source sensitivity compared to slow scan speed of 10 arcsec/sec, 
while not yet causing point spread function (PSF) degradation 
which is present in fast scans at 60 arcsec/sec because of data 
sampling and detector time constants. 

Much of the science from deep far-infrared surveys strongly 
benefits from the availability of ancillary multiwavelength data 
that historically have been obtained in fields of a certain shape 
and orientation. We have matched our blank field data to such 
constraints (e.g. the shape/orientation of the mid-IR GOODS 
fields) by typically defining scan layouts in 'sky' coordinates 
rather than 'array' coordinates that rotate with epoch of observa- 
tion. This is also facilitating the combination of data from differ- 
ent epochs. Exceptions to this approach are the COSMOS field 
(discussed below) and the small lensing cluster fields for which 
there is no preferred orientation. 

Reaching the desired depths for our fields requires multiple 
passages over a given point in the sky. Rather than simply re- 
peating a single scan we use this to improve the redundancy of 
the data. Specifically, the observation setup for PEP includes 

- Scans in both nominal and orthogonal directions. While the 
highpass-filtered data reduction discussed in Sect. [4] does 
not require this, such crosslinking is essential for alterna- 
tive reductions with full inversion algorithms such as the 
MadMap implementation in the Herschel Common Science 
System (HCSSfl software environment. Typically, we obtain 
the same number of astronomical observing requests (AORs) 
in nominal and orthogonal scan direction, often concatenated 
in pairs. For the very elongated EGS and z~l clusters we 
overweight maps scanning along the long axis, in order to 
reduce overhead losses that are caused by the time spent in 
scan turnaround loops. 

1 HCSS is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia 
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Table 2. AOR parameters used for the PEP fields 
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Notes. For both the nominal and the orthogonal scan direction, the table lists the scan leg length, the cross-scan separation step, the number of 
scan legs, the scan angle and its reference system (Sky or Array), the number of scan repetitions within an AOR and the number of independent 
AORs obtained with those settings. 'Horn Sq' indicates that the cross-scan separation and number of scan legs are defined automatically by the 
AOT logic to produce a square map with 'homogeneous' coverage. 



- Small cross-scan separations, specifically the size of one of 03 00 
the eight PACS 'blue' detector array matrices (~50") or frac- 
tions 1/2 or 1/2.5 of it. Simple models demonstrate that for 
such a scan pattern homogeneous coverage maps are pro- 45 
duced already from a single detector matrix, for any relative 
orientation of the scan direction and the inter-matrix gaps. 
By definition, this pattern also averages out sensitivity varia- 3c 
tions between detector matrices over most of the final map. If 
this redundant mapping scheme led to short execution times 
of a single scanmap over the field but deeper observations ° 15 
were needed, the scanmap was repeated within an AOR to | 
reach a total execution time between one and a few hours. a 

- Often, many AOR pairs with a plausible execution time that 02 o 00 , 
is not exceeding a few hours each, are still needed to achieve 
the required depth. Then, AOR positions may be dithered 
by a fraction of the cross-scan separation to further improve 
spatial redundancy, and the corner of the map where the scan 
is started may be varied. 



The two square degree COSMOS field is a special case 
where the described observing strategy would lead to extremely 
long individual AORs. During one of the about two month long 
Herschel visibility periods of COSMOS, the PACS arrays are 
always similarly oriented on this region of the sky, with a posi- 
tion angle of the long axis about 20±10 °. This permits an AOR 
setup in the efficient 'homogeneous coverage' mode in array co- 
ordinates, keeping the individual AOR length below 5 hours but 
staying well matched to the roughly square nonrotated orienta- 
tion of many COSMOS ancillary data sets. This is illustrated by 
the actual coverage map obtained with PACS (Fig. [2]). 

Table|2]provides the key parameters of the actual PACS AOR 
implementations for each field. No significant source variability 
is expected for the dust-dominated emission of almost all de- 
tected sources. For this reason no timing contraints needed to 
be applied in the scheduling. For practical reasons scheduling of 
all AORs of a field during a visibility period was aimed for, and 
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Fig. 2. Actual 160 /itn coverage map of the COSMOS field. 
Spots of reduced coverage near the top left of the map are due to 
dropping of data because of 'speed bumps' (see text). 

was typically but not always achieved. For fields near the plane 
of the ecliptic (COSMOS), asteroid passages may introduce an- 
other time dependent factor, as clearly demonstrated in mid- 
infrared detections during Spitzer observations of the COSMOS 
field (Sanders et al. 2007 ). The contrast between galaxies and as- 
teroids is more favourable in the far-infrared. Still, bright aster- 
oids would be detectable in individual maps if present but were 
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Fig. 3. Offsets of GOODS-N 100 /mi sources detected at > 5<x 
in a 'blind' extraction, compared to the corresponding 24 jjm 
positions of the MIPS catalog used by Magnelli et al. (|2009). 
A global offset has already been subtracted when deriving the 
PACS map. 

not identified when differencing our individual COSMOS maps 
from a coaddition. 

SPIRE maps for most of the PEP fields are obtained in co- 
ordinated observations by the HerMES key program (Oliver et 
al. 1201 11 1. For the two z~l clusters we implemented within PEP 
simple 10'xlO' 'large' SPIRE scanmaps in nominal scan speed, 
spatially dithering between five concatenated independent repe- 
titions. 



4. Data Analysis 

4. 1 . Reduction of scanmap data and map creation 

For scanning instruments with detectors that have a significant 
1/f low frequency noise component, map creation usually fol- 
lows one of two alternative routes. One is using full 'inversion' 
algorithms as widely applied by the cosmic microwave back- 
ground community and the other uses highpass filtering of the 
detector timelines and subsequent direct projection, frequently 
used for Spitzer MIPS 70 or 160 /an reductions. An algorithm 
of the first 'inversion' type is available in the HCSS Herschel 
data processing in the form of an implementation and adaption 
to Herschel of a version of the MadMap code (Cantalupo et al. 
120101) . The alternative option that we adopt and describe in more 
detail below is using highpass filtering of the detector timelines 
and a direct 'naive' mapmaking. This choice is made because for 
our particular case of deep field observations, MadMap presently 
does not reach the same point source sensitivity, and the preser- 
vation of diffuse emission is not important for our science case. 
As noted, the cross-linked design of the PEP observations how- 
ever does permit the future application of such inversion codes. 

Our reduction first proceeds on a per AOR level and is 
based on scanmap scripts for the PACS photometer pipeline 
(Wieprecht et al. 1200 9) in HCSS, with parameter settings and 



additions optimized for our science case. After retrieving PACS 
data and satellite pointing information we apply the first reduc- 
tion steps to the time-ordered PACS data frames, identifying 
functional blocks in the data, flagging bad pixels, flagging any 
saturated data, converting detector signals from digital units to 
volts and the chopper position from digital units to physical an- 
gle. 

After adding to the time-ordered data frames the instanta- 
neous pointing obtained from the Herschel pointing product, we 
apply 'recentering' corrections. These are derived by comparing 
PACS maps obtained in a separate first processing of partial dat- 
ablocks to deep 24 //m catalogs with accurate astrometry. Deep 
radio catalogs can also be used successfully. The corrections are 
derived by stacking at the positions of 24 /mi sources a PACS 
100 /mi or 70 //m map that is obtained from typically about 15 
minutes of data, and measuring the offset of the stacked PACS 
detection. The measured position offset is then corrected for this 
block of data in the original frames. The blocks used were typi- 
cally restricted to one scan direction (all 'odd' or 'even' scanlegs 
of a map repetition, or subsets for the large COSMOS maps). By 
this procedure we correct for (1) the global offset pointing error 
of Herschel which is of order 2" RMS and thus nonnegligible 
compared to PACS beamsizes, (2) differential pointing errors be- 
tween different AORs that would lead to PSF smearing, in partic- 
ular if AORs have been observed spread over long periods with 
different spacecraft orientations, (3) residual timing offsets be- 
tween PACS data and pointing information that manifest them- 
selves in small offsets between odd and even numbered scan- 
lines and (4) pointing drifts over >0.5hour timescales within an 
AOR. Overall this leads to subarcsecond astrometry and reduces 
unnecessary PSF smearing by pointing effects. As a verification 
of this approach, Figure [3] compares the positions from a blind 
100 //m catalog based on the final GOODS-N map with 24 //m 
positions. 

Short glitches in the detector timelines caused by ionizing 
particle hits are flagged and interpolated with an HCSS imple- 
mentation of the multi-resolution median transform, developed 
by Starck & Murtagh d 19981 ) to detect faint sources in ISOCAM 
data. The different signatures of real sources and glitches in 
the pixel timeline at medium scan speed are discriminated by a 
multi-scale transform. This procedure is known to produce false 
detections on bright point sources. This was no problem for our 
deep field data except for a few bright COSMOS sources for 
which we locally reduced the glitch detection sensitivity. 

The '1/f noise of the PACS photometers, in fact roughly 
oc /~ - 5 over the relevant frequencies, is removed by subtracting 
from each timeline the timeline filtered by a running box median 
filter of radius 15 samples (30") at 70 fim or 100 /mi and 26 sam- 
ples (52") at 160 fim. A mask is used to exclude sources from 
the median derivation. The mask is created by thresholding a S/N 
map produced from a smoothed coadded map that is including 
all AORs. Tests were done adding simulated sources to the time- 
lines before masking and before highpass filtering. These tests 
were done on the real timelines of the full Lockman Hole AOR 
set for the red filter, thus implicitly including all real sky struc- 
tures. They indicate that the filtering modifies the fluxes by less 
than 5% for masked point sources and ~16% for very faint un- 
masked point sources. These results apply to 'extragalactic deep 
field' skies that are composed of many point sources and to our 
specific highpass filter setting. Extended sources need different 
filtering radii or reduction methods. 

Flatfielding and flux calibration is done via the standard 
PACS pipeline calibration files. Because of the excellent in-flight 
stability of the photometer response, no use is made of the ob- 
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Fig. 4. Correlation map extracted from GOODS-N 100 11m data, 
reflecting correlations induced by map projection and - in scan 
direction - by residual low frequency noise and highpass filter- 
ing. Correlations across the GOODS-N map are highly uniform 
and depend only on the relative Ax, Ay pixel positions, repre- 
sented here with respect to the central pixel. 



servations of the calibration sources internal to PACS that are 
obtained within each AOR. PEP science demonstration phase 
data were reduced using earlier versions of the responsivity cali- 
bration file and corrections derived on the stars yDra, aTau, and 
aCMA; they are within 5% of the currently valid Version 5 of the 
reponsivity calibration file that has been validated over a wider 
set of flux calibration stars and asteroids. 

Before creating maps from the timelines we discard various 
unsuitable parts of the data: 

1. Observations of the internal calibration sources and a short 
subsequent interval during which the signal re-stabilizes. 
SDP observations had an erroneously large frequency of 
such calibration blocks. 

2. Data taken during the turnaround loops between scans. The 
nominal pointing accuracy is not guaranteed during this 
phase, and highpass filtering effects on source flux are more 
severe if the telescope slows down strongly or even stops 
during a subtantial part of the filtering window. 

3. Before a reduction of the star tracker operating temperature 
in Herschel operational day 320, 'speed bumps' in the satel- 
lite movement occured, when stars crossed irregularly be- 
having star tracker pixels. This led to an unkown mispoint- 
ing. These events can be identified by comparing the two 
angular velocities obtained from (1) differentiating the po- 
sitions on sky reported in subsequent entries of the point- 
ing product and (2) the direct angular rate information in the 
pointing product. 

4. Blue PACS bolometer data occasionally show fringes in 
the signal due to magnetic field interference that is heavily 
aliased into the timelines and final maps. Maps created sep- 
arately for each individual scan line were visually inspected 
and scans clearly affected by fringing excluded from final 
mapping. 

Maps are obtained from the timelines for each AOR via the 
HCSS 'photProject' projection algorithm, which is equivalent to 
a simplified version of the 'drizzle' method (Fruchter & Hook 
120021 ). Given the high data redundancy in the deep fields, PSF 
widths and noise correlation in the final map can be reduced 
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Fig. 6. Results of artifical source addition experiments for the 
GOODS-N 160 fim extractions. Top: prior extraction, Bottom: 
blind extraction. The red lines represent the average photometric 
accuracy, blue lines set the standard deviation observed in each 
flux bin (after 3<x clipping). The systematic boosting in the blind 
extraction is corrected in the final blind catalog. Given that, in- 
put and output fluxes are consistent with each other within a few 
percent, testifying that PACS fluxes are reliably extracted and 
retrieved by the adopted source extraction methods. The orange 
histogram represents the detection rate (or completeness) com- 
puted on the artificial injected sources, while black histograms 
denote cp the fraction of spurious sources as a function of out- 
put flux. 



choosing smaller projection drops than the physical PACS pixel 
size. Drop sizes between 1/8 and 1/4 of the physical PACS pixel 
size were used, depending on redundancy of the individual maps 
for each field. Weights of the different detectors in the projec- 
tion consider the inverse variance derived from the noise in the 
dataset itself. 

Maps from each AOR were coadded into final maps, weight- 
ing the indivudal maps by the effective exposure of each pixel. 
The final error map was computed as the standard deviation of 
the weighted mean. PSF fitting using the methods described in 
Sect. 14.21 assumes errors that are uncorrelated between neigh- 
bouring pixels. In practice, correlations exist due to projection 
and due to to the correlations that are caused by residual 1/f 
noise in the filtered timelines, in particular along the scan di- 
rection. We have verified that, because of the high redundancy 
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Fig. 5. (Left) PACS 100 pm to MIPS 24 pm flux density ratio as a function of the PACS 100 /jm prior catalog flux density. Shaded grey regions 
show the space density distribution of galaxies observed in the GOODS-S field. In that field, the MIPS 24 pm and PACS 100 pm catalogs reach at 
3-cr limit of 20 p]y and 1.2 mjy, respectively. Red dots and orange triangles present the evolution with redshift of the PACS-to-MIPS flux density 
ratio of ultra-luminous (ULIRGs) and luminous (LIRGs) infrared galaxies, respectively, predicted using the Chary & Elbaz (2001 ) library. Each 
symbol corresponds to a given redshift, in intervals of Az = 0.5. On each track, we indicate on one point its corresponding redshift as well as as the 
path followed for increasing redshifts. The dashed black line represents the limit of the parameter space reachable using the MIPS 24 pm catalog 
available in this field; all the PACS 100 pm sources below this line would have a MIPS 24 pm counterpart in this catalog. (Right) PACS 160 pm to 
MIPS 24 /jm flux density ratio as function of the PACS 160 /jm flux density. Symbols and lines are the same as in the left panel. 
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Fig. 7. Noise estimate from randomly placing apertures on the 
GOODS-N 160 pm residual map. 



of the data, these correlations are close to uniform across the fi- 
nal map, with less than 2% variation on the correction factor that 
is derived below. Thus we derived from PSF shape and correla- 
tion information a mean correlation correction factor which was 
then accounted for in the errors on the extracted fluxes. A corre- 
lation map is constructed, starting by collecting series of paired 
pixel values with same relative pixel coordinate offsets i,j. We 
take values Af defined as the deviation in flux of a pixel in an 
individual AOR map from the corresponding pixel in the final 
map. This is a deviation from the mean with expectation value 
0. Such values are taken for a large number of pairs in different 
positions in each AOR map and from different AORs. For each 
pixel-pairs series Afj.2, which corresponds to a specific i,j pixel 
offset, a correlation coefficient is calculated: 



P(i, j) = 



EA/iA/2 



V2(A/i) 2 V2(A/ 2 ) 2 



(i) 



The correlation coefficients are stored as a map and written to the 
i,j position relative to the central pixel. Figure [4] shows an exam- 
ple correlation map. Normal error propagation for a weighted 
sum g{x\..x„) = Y!l=i a k ' x ki where the error are cru and correla- 
tion coefficients p(k, I) is: 



^ a k o- k ■ a,cri ■ p(k, I) 



(2) 



k,l=\ 



Given a PSF stamp with pixel values P k and a correlation be- 
tween every two pixels p(i, j) known from their relative position, 
the correlation correction factor to the derived errors is: 



YjPiPjp(iJ) 



(3) 



/ is the ratio of the propagated error without correlations and the 
error calculated with the correlation terms. This assumes a near 
uniform error map on a scale of a PSF. For the typical pixel sizes 
(2" and 3" at 100 pm and 160 pm) and projection parameters 
used in the PEP SDP reductions, the corrections are about / « 
1 .4 for the 70 pm and 100 pm maps and / « 1 .6 for the 160 pm 
maps. 

Before availability of version 6 of the Arraylnstrument cal- 
file that is containing the spatial transformation from PACS fo- 
cal plane to sky, 160 pm data showed as 1" spatial offset from 
70 pm or 100 pm data. We corrected for this ad hoc in the 160 pm 
map fits headers, using offsets derived from a comparison of pre- 
liminary catalogs in the two bands. 

4.2. Catalog creation 

We extract source catalogs using point source fitting routines 
outside HCSS. We have used both blind extraction via the 
Starfinder PSF-fitting code (Diolaiti et al. 2000) and a guided 
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extraction using 24 fim source priors, following the method de- 
scribed in Magnelli et al. (2009). We fitted with point spread 
functions extracted from the maps. Since these observed PSFs 
are limited in radius, we used for aperture corrections point 
spread functions obtained on Vestal In order to match the obser- 
vations, these were rotated to match the satellite position angle 
for each observation of a field, coadded and slightly convolved 
with a gaussian to match the actual FWHM of the combined 
map. For our reduction methods and projection into 2" map pix- 
els at 70 and 100 fim and 3" at 160 fim, we have PSF FWHM 
for the 70, 100, and 160 fim maps of 6.46", 7.39+0.10", and 
1 1.29+0. 1". The FWHM values are from gaussian fits to the 
core of the observed PSF. For the 100 and 160 fim widths we 
quote the error of the mean of measurements from five different 
fields. 

All PEP blank fields benefit from extensive multi- 
wavelength coverage that is allowing guided extraction based 
on source positions at shorter wavelengths, where the depth and 
resolution of the observations are higher. This approach resolves 
most of the blending issues encountered in dense fields and al- 
lows straightforward multi-wavelength association (Magnelli et 
al. 120091 l20m Roseboom et al. IMTOb . However, to use this 
powerful method, prior source catalogs have to contain all the 
sources in the PACS images. Deep MIPS 24 fim observations, 
available for all our blank fields, should fulfill this criterion 
since they have higher resolution and are deeper than our cur- 
rent PACS observations. Moreover, since the 24 fim emission is 
strongly correlated with the far-infrared emission, those catalogs 
will not contain a large excess of sources without far-infrared 
counterparts. This largely avoids deblending far-infrared sources 
into several unrealistic counterparts, as could happen when using 
an optical prior catalog with very high source density. Figure [5] 
illustrates the validity of the MIPS 24 /mi observations as PACS 
prior source positions. 

At 100 /urn, models and observations predict typical PACS- 
to-MIPS flux ratios in the range 5-50. A MIPS 24 fim catalog 50 
times deeper than the PACS observations is thus suited to per- 
form guided source extraction. This is illustrated, for GOODS-S, 
in the left panel of Figure [5] where the observed PACS popula- 
tion lies well-below the boundary of the parameter space reach- 
able by the MIPS 24 /mi catalog. In all our blank fields, deep 
MIPS 24 /zm catalogs fulfill this criterion (i.e., GOODS-S/N, 
COSMOS, LH, ECDFS, EGS. . . ). 

At 160 //m, models and observations have typical PACS-to- 
MIPS flux ratios in the range 15-150. In all but one field, the deep 
MIPS 24 fim observations are at least 150 times deeper than our 
PACS observations; thus, they can be used as prior source posi- 
tions. In GOODS-S, the deepest MIPS 24 fim observations are 
only 100 times deeper than our PACS observations. This limita- 
tion can be observed, in the right panel of Figure |5j as a slight 
truncation, at faint 160 /mi flux density, of the high-end of the 
dispersion of the PACS-to-MIPS flux ratio. This truncation will 
introduce, at faint flux density, incompleteness in our prior cata- 
log. However, we also observe that even at this faint 160 fim flux 
density, the bulk of the population has a PACS-to-MIPS flux ra- 
tio of ~ 60. The incompleteness introduced by the lack of MIPS 
24 fim priors should thus be low or at least lower than the in- 
completeness introduced by source extraction methods at such 
low S/N. This was checked by comparing the GOODS-S 160 
fim catalog obtained using blind source extraction with that ob- 
tained using guided source extraction: we find no significant dif- 
ference, at faint 160 fim flux densities, in the number of sources 
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in those two catalogs. In line with these finding, Magdis et al. 
(in prep) find very low fractions of 24/zm undetected sources in 
very deep data from the GOODS-Herschel key program, < 2% 
for ratios of detection limits S(100)/S(24)>43 and < 1% for 
S(160)/S(24)>130. 

Therefore, for all fields with deep MIPS 24 fim observations, 
we also extract source catalogs with a PSF-fitting method using 
24 fim source positions as priors and following the method de- 
scribed in Magnelli et al. (2009). We used the same PSFs and 
aperture corrections as for the blind source extraction. Blind and 
prior catalogs were compared to verify the consistency between 
those two methods. 

Completeness, fraction of spurious sources and flux relia- 
bility were estimated by running Monte Carlo simulations. Up 
to 10000 artificial sources were added to PACS science maps, 
and then extracted with the same techniques and configurations 
adopted for real source extraction. In order to avoid crowding, 
many such frames were created, each including a limited num- 
ber of artificial sources. The number of frames and the number 
of sources added in each one depend on the size of the field un- 
der analysis and range between 20 and 500 sources (GOODS 
fields or COSMOS) per frame, repeated up to reaching the total 
of 10000. These synthetic sources cover a large range in flux, 
extending down to 0.5cr (cr being the measured rms noise in the 
PACS maps). The flux distribution follows the detected num- 
ber counts, extrapolated to fainter level by means of the most 
successful fitting backward evolutionary model predictions (see 
Berta et al. 1201 11 1. Sources are modelled using the Vesta PSF, 
manipulated as described above. 

Figure [6] shows an example of results in the GOODS-N 
field at 160 fim. Completeness is defined here as the fraction 
of sources that have been detected with a photometric accuracy 
of at least 50% (Papovich et al. 120041 ). Spurious sources are de- 
fined as those extracted above 3<x with an input flux lower than 
3cr(Image). The systematic flux boosting in the blind extraction 
is corrected in the final blind catalog on the basis of these simu- 
lations. 

Noise was estimated by extracting fluxes through 10000 
apertures randomly positioned on residual maps. FigureQshows 
the distribution of the extracted fluxes, peaking around zero, 
as expected for a well subtracted background, and showing a 
Gaussian profile. 

5. Science demonstration phase data 

Figures [8] and [9] show the 100 fim and 160 fim maps of 
the GOODS-N and Abell 2218 fields as obtained during the 
Herschel science demonstration phase. A conservative thresh- 
old of 90% completeness is reached for the blind catalogs near 
7 mJy and 15 mJy for 100 fim and 160 fim in the main parts 
of both GOODS-N and A2218. Above that level, the blind cata- 
logs contain 153 and 126 sources for GOODS-N and 49 and 47 
for Abell 2218, respectively. At the time of these observations, 
the Herschel scan maps were still exhibiting larger turnaround 
overheads than implemented later, and the PACS scanmap AOR 
had an unnecessarily large frequency of internal calibrations. 
Both these factors do not significantly affect the maps which are 
based on highpass-filtered reductions. With the exception of the 
overheads, these observations are representative for the results 
achievable in later mission periods during the observing times 
listed in TableQ] which reflect the later reduced overheads. 

Figure [TOl compares 100 fim fluxes for the GOODS-N field 
between the blind (starfinder) extraction and extraction based on 
24 fim priors. The blind detections have here been associated a 
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Fig. 9. Science demonstration phase maps of the Abell 2218 lensing cluster field. Left: 100 fim. Right: 160 fim 



posteriori to 24 fim sources, using a modification of the maxi- 
mum likelihood method of Ciliegi et al. d20011 l. The agreement 
is satisfactory with no systematic flux differences. Deviations oc- 
cur at low fluxes where either catalog is incomplete and for few 
outliers where algorithms disagree in splitting a peak into two 
sources compared to one. 

Figure Q~T| shows a 'color-magnitude' diagram for 160 fim 
detected sources in GOODS-N. As expected, there is a tendency 
for higher redshift sources to be redder in the 160/100 fim flux 



ratio, with considerable scatter due to measurement error and 
variation among the population. A detailed discussion of col- 
ors and SEDs is outside the scope of this work, see also Elbaz 
et al. d20T0b . Mag nelli e t al. d20T0T l. Hwang et al. d20T0b and 
Dannerbauer et al. (2010) for first results on GOODS-N spectral 
energy distributions. 

Figure[l2]illustrates the potential of PEP to reach normal star 
forming galaxies up to redshifts z~2. The GOODS-N science 
demonstration phase data reach at z~2 infrared luminosities of 
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Fig. 10. Comparison of blindly extracted and prior extracted 
100 fim fluxes for the GOODS-N field. 
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Fig. 12. Red: Total infrared luminosity of GOODS-N sources de- 
tected above 3<x in the PEP SDP data. Blue: Sources from the 
GOODS-S and COSMOS data illustrate the potential of the full 
PEP survey. 



As an example of the PEP/PACS performance for deep sur- 
veys, the > 90% completeness blind catalogs of the GOODSN 
and Abell 2218 SDP data will be released on the PEP websiteB 
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10 12 L G ( s,;ar formation rates about 100 solar masses per year). 
The GOODS-N SDP observation thus reach the star formation 
'main sequence' for massive galaxies at that redshift (Daddi et 
al. I2007ab . The GOODS-S data go deeper by a factor >2, and 
the COSMOS data provide the precious statistics at higher lumi- 
nosities. 



6. Overview of first science results of the PEP 
survey 

6. 1 . The cosmic infrared background 

One of the most imminent tasks of a large far-infrared space tele- 
scope is to resolve a large fraction of the cosmic infrared back- 
ground into its constituent galaxies. The CIB peaks at roughly 
150 fim (e.g. Dole et al. l2U0"6l making Herschel/PACS excel- 
lently suited to characterise the sources of its bulk energy out- 
put. Berta et al. d20T0l l used SDP observations of GOODS-N to 
derive 100 and 160 fim number counts down to 3.0 and 5.7 mjy, 
respectively. Altieri et al. ( 1201 Oi l used strong lensing by the mas- 
sive cluster Abell 2218 to push this limit down to 1 and 2 mjy, 
deeper by a factor ~3 for both wavelengths. Berta et al. (1201 11 1 
use a complement of PEP blank field observations including the 
deep GOODS-S observations to reach a similar depth as Altieri 
et al. but with the larger statistics provided by the blank field ob- 
servations. They also add 70 fim counts overall reaching 1.1, 1.2, 
and 2.4mJy in the three PACS bands. To t his dep th, 58% (74%) 
of the COBE CIB as quoted by Dole et al. (12006b is resolved into 
individually detected sources at 100 (160 fim). These fractions 
reach 65% (89%) when including a P(D) analysis. 

It is important to note that the direct COBE/DIRBE mea- 
surements have considerable uncertainties due to the difficulty 
of foreground subtraction. Similar to other wavelengths such as 
the <10keV X-rays (e.g. Brandt & Hasinger 2006) we are reach- 
ing the point where the lower limits to the cosmic background 
that are provided by the integral of the resolved Herschel mea- 
surements are more constraining than the direct measurements. 
Quoting 'resolved fractions' in reference to COBE hence starts 
to become problematic. 

Because of the excellent multiwavelength coverage of the 
PEP fields, spectroscopic or photometric redshifts can be as- 
signed to the detected sources, allowing us to determine the red- 
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shifts at which the CIB originates. To the depth reached by Berta 
et al. d20TTD , half of the resolved CIB originates at z>0.58, 0.67, 
0.73 for wavelengths of 70, 100, and 160 /mi, respectively. These 
redshifts are mild lower limits because they exclude faint sources 
that are not individually detected by Herschel, and naturally in- 
crease with wavelength, reflecting the dominant contributions by 
the far-infrared SED peak of increasingly distant objects. 

Comparing the PEP counts and redshift distributions to back- 
ward evolutionary models tuned on pre-Herschel data, Berta et 
al. ( 1201 11 1 find for several models reasonable agreement with the 
total counts but more significant mismatches to the observed red- 
shift distributions. Clearly, models have to be tuned in order to 
provide a satisfactory representation of the most recent data in- 
cluding Herschel. 

The PEP GOODS-N SDP data have been used to derive 
first direct far-infrared based luminosity functions up to z~2-3 
(Gruppioni et al. 1201 Oi l. Strong evolution of the comoving in- 
frared luminosity density (proportional to the star formation rate 
density) is found, increasing with redshift as (1 + 2) 3 - 8±0 - 3 up to 
z~l. Global classification of the SEDs assigns to most objects 
either a starburst-like SED or an SED that is suggesting a mod- 
est AGN contribution. 

6.2. 'Calorimetric' far-infrared star formation rates and star 
formation indicators 

Star formation rates are one of the key measurables in galaxy 
evolution studies. 'Calorimetric' rest frame far-infrared mea- 
surements are the method of choice for massive and in partic- 
ular dusty galaxies but have been out of reach for typical high- 
z galaxies in the pre-Herschel situation. Typically, star forma- 
tion rates were measured from the rest frame ultraviolet, mid- 
infrared, submm/radio, or a combination of those methods. All 
of them require certain assumptions. Extracting a star formation 
rate from the rest frame ultraviolet continuum involves breaking 
the degeneracies between the star formation history and obscu- 
ration, and involves assumptions about the dust extinction law 
and/or geometry of the obscuring dust. Extrapolation from the 
mid-infrared or submm/radio requires the adoption of SED tem- 
plates or assumptions that were derived locally but insufficiently 
tested at high redshift. 

Nordon et al. (1201 Oi l focussed on massive z~2 normal star 
forming galaxies that are currently the subject of intense study 
towards the role of secular and merging processes in their evolu- 
tion (e.g. Forster Schreiber et al. 2009 ). They found ultraviolet- 
based star formation rates reasonably confirmed. At the rest 
frame optical bright end of these bright star forming galaxies, 
UV SFRs are overpredicting the calorimetric measurement by a 
factor 2 only. In contrast, a significant overprediction by a factor 
~4- 7.5 was found when extrapolating from the 24 /mi flux as- 
suming that the local Universe Chary & Elbaz (2001 ) templates 
apply for the given mid-IR luminosity. This is in line with the 
encompassing SED studies of Elbaz et al. (2010) and Hwang et 
al. (|2010) which, combining PEP and HerMES as well as lo- 
cal Akari data, found that this overprediction by 24 yum extrap- 
olation sets in at z~1.5, and that the FIR SED temperatures of 
high z galaxies are modestly colder than their local equivalents 
of similar total infrared luminosity. 

There are two ways of looking at these results. Comparing 
at same total IR luminosity local galaxies with z~2 galaxies, the 
z~2 'overprediction' when extrapolating from the mid-IR could 
either be due to an enhanced ratio of the PAH complex to the 
FIR peak, or due to a boosting of the mid-infrared emission via 
the continuum of a possibly obscured AGN (see also Papovich et 



al. 120071 and Daddi et al. l2007bl ). While most likely both factors 
are at work to some level, the setting in of the 24 /mi overpredic- 
tion at z~1.5 as well as Spitzer spectroscopy of mid-IR galaxies 
at these redshifts (Murphy et al. 2009 Fadda et al. 120101 ) ar- 
gues for a dominant role of stronger PAH. This is established by 
Nordon et al. 1201 II who show that z~2 GOODS-S galaxies with 
large 24 /mi/FIR ratios indeed have PAH dominated spectra in 
the Fadda et al. (2010) sample. 

It is important to put this into the perspective of other prop- 
erties of Lir ~ 10 12 L o local and z~2 galaxies. Locally, such 
an object is a classical interacting or merging 'ultraluminous 
infrared galaxy' (ULIRG, Sanders & Mirabel 1996) with star 
formation concentrated in a few hundred parsec sized region. 
Such star formation rates are well above the local 'main se- 
quence' (Brinchmann et al. 2004). At z~2, objects with the same 
star formation rates can be on the main sequence (Daddi et al. 
2007a ) and are often not interacting/merging but massive, turbu- 
lent disks with star formation spread out over several kpc (e.g. 
Genze l et al. 120081 Shapiro et al. 120081 Forster Schreiber et al. 
2009). It is hence not surprising that SED templates calibrated 
on local compact mergers fail to reproduce the more extended 
star formation at z~2. The locally defined connotations of the 
'ULIRG' term, beyond its basic definition by IR luminosity, ob- 
viously do not apply automatically at high redshift. 

Rodighiero et al. (fJOlO) use the direct Herschel far-infrared 
star formation rates in the PEP GOODS-N field to characterize 
the evolution of the Specific Star Formation Rate (SSFR)/mass 
relation up to z~2. They find a steepening from an almost flat 
local relation to a slope of -0.5 at z~2. 

Concerning the most highly star forming high redshift ob- 
jects, Magnelli et al. d20T0b combine the GOODS-N and Abell 
2218 PACS data with submm surveys. They use PACS to sample 
the Wien side of the far-infrared SED of submm galaxies and 
optically faint radio galaxies with accurately known redshifts. 
The directly measured dust temperatures and infrared luminosi- 
ties are in good agreement with estimates that are based on radio 
and submm data and are adopting the local universe radio/far- 
infrared correlation (see also the dedicated PEP/HerMES test of 
the high-z radio/far-infrared correlation by Ivison et al. ( 120101 ). 
This confirmation of huge O(lOOOM yr _1 ) star formation rates 
for SMGs is in support of a predominantly merger nature of 
SMGs, since such star formation rates are very hard to sustain 
with secular processes (e.g. Dave et al. 120 101 ). 

6.3. The role of environment 

PEP data cover a wide range of environments. As first steps, 
Magliocchetti et al. (2011) derive correlation functions and co- 
moving correlation lengths at z~l and 2 from the GOODS-S 
data which also provide evidence for most infrared bright z~2 
galaxies in this field residing in a filamentary structure. At z~l, 
Popesso et al. d201 II ) observe a reversal of the local star forma- 
tion rate - density relation which is linked to the presence of 
AGN hosts that are favoring high stellar masses, dense regions 
and high star formation rates. 

6.4. Studies of individual galaxy populations 

Dannerbauer et al. (120101) use the < 10" beam PEP maps to 
verify the identifications of (sub)millimeter galaxies and test 
the potential of PACS colors and mid-infrared to radio photo- 
metric redshifts for these objects. Santini et al. ( 120101 ) use the 
Herschel points to break the degeneracy between dust temper- 
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ature and dust mass that is inherent to submm-only data, and 
derive high dust masses for a sample of submillimeter galaxies. 
Comparison of these high dust masses with the relatively low 
gas phase metallicities either implies incompletely understood 
dust properties or a layered structure which is combining low 
metallicity visible outer regions with a highly obscured interior. 

Magdis et al. (J2010) use a stacked detection of z~3 Lyman 
break galaxies to constrain their far-infrared properties. Similar 
to the SMG/OFRG study of Magnelli et al. (l2l)T0l ). this work 
highlights the selection effects in the luminosity / dust tem- 
perature plane that are imposed by groundbased submm detec- 
tion. Bongiovanni et al. (120101) study Lyman a emitters in the 
GOODS-N field. PACS detections of part of these galaxies are 
evidence for overlap with the dusty high redshift galaxy popula- 
tion. 

6.5. AGN-host coevolution: Two modes? 

Within the complex task of disentangling the coevolution of 
AGN and their host galaxies, and the role of AGN feedback, rest 
frame far-infrared observations provide a unique opportunity to 
determine host star formation rates. This rests on the assump- 
tion that the host dominates over the AGN in the rest frame far- 
infrared. This assumption is motivated by determinations of the 
intrinsic SED of the AGN proper, which is found to drop towards 
the far-infrared (e.g. Netzer et al. 120071 Mullaney et al. 1201 11 1. 
Compared to previous attempts from the submm (e.g. Lutz et 
al. 1201 Oi l and Spitzer (e.g. Mullaney et al. 1201 Oi l. He rschel per- 
formance provides a big step forward. Shao et al. ( 20 1 ()|l use 
the PEP GOODS-N observations and the Chandra Deep Field 
North to map out the host star formation of AGN of different 
redshift and AGN luminosity. The host far-infrared luminosity 
of AGN with L2-iokeV ~ 10 43 erg s -1 increases with redshift by 
an order of magnitude from z=0 to z~l, similar to the increase 
with redshift in star formation rate of inactive massive galaxies. 
In contrast, there is little dependence of far-infrared luminosity 
on AGN luminosity, for La-iokeV S 10 44 erg s _1 AGN at z>l. In 
conjunction with properties of local and luminous high-z AGN, 
this suggests an interplay between two paths of AGN/host coevo- 
lution. A correlation of AGN luminosity and host star formation 
for luminous AGN reflects an evolutionary connection, likely via 
merging. For lower AGN luminosities, star formation is similar 
to that in non-active massive galaxies and shows little depen- 
dence on AGN luminosity. The level of this secular, non-merger 
driven star formation increasingly dominates over the correlation 
at increasing redshift. 

7. Conclusions 

Deep Herschel far-infrared surveys are a powerful new tool for 
the study of galaxy evolution and for unraveling the constituents 
of the cosmic infrared background. We have described here the 
motivation, field selection, observing strategy and data analy- 
sis for the PEP guaranteed time survey. Science demonstration 
phase data of GOODS-N and Abell 2218 are discussed to illus- 
trate the performance of Herschel-PACS for this science. The 
wide range of initial science results from the PEP data is briefly 
reviewed. 
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